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ABSTRACT 

We investigate effects of baryon dissipation on the dark matter virial scaling relation between total mass and 
velocity dispersion and the velocity bias of galaxies in groups and clusters using self-consistent cosmological 
simulations. We show that the baryon dissipation increases the velocity dispersion of dark matter within the 
virial radius by « 5 - 10%. The effect is mainly driven by the change in density and gravitational potential in 
inner regions of cluster, and it is larger in lower mass systems where gas cooling and star formation are more 
efficient. We also show that the galaxy velocity bias depends on how galaxies are selected. Galaxies selected 
based on their stellar mass exhibit no velocity bias, while galaxies selected based on their total mass show 
positive bias of « 10%, consistent with previous results based on collisionless dark matter-only simulations. 
We further find that observational estimates of galaxy velocity dispersion are unbiased with respect to the ve- 
locity dispersion of dark matter, provided galaxies are selected using their stellar masses and and their velocity 
dispersions are computed with more than twenty most massive galaxies. Velocity dispersions estimated with 
fewer galaxies, on the other hand, can lead to significant underestimate of dynamical masses. Results pre- 
sented in this paper should be useful in interpretating high-redshift groups and clusters as well as cosmological 
constraints derived from upcoming optical cluster surveys. 

Subject headings: cosmology: theory-clusters: formation- methods: numerical 



1. INTRODUCTION 

Galaxy clusters are unique laboratories for astrophysics and 
powerful probes in cosmology. As the largest and most mas- 
sive gravitationally-bound structures in the universe, they pro- 
vide well-defined environments for studies of various astro- 
physical processes: e.g., plasma interactions in the intraclus- 
ter medium, galaxy interactions and evolution, supermassive 
black hole feedback processes, etc. In the cosmological con- 
text, galaxy clusters are good tracers of the cosmological evo- 
lution, as thei r ov erall dyn amics are largely driven by gravity 
(Kaiser 1986, see Voit 2005 for a review). Different ways in 
which clusters can be used for cosmology rely on understand- 
ing of how their observable properties are related to the total 
virialized mass. One of the main challenges of cluster cos- 
mology and astrophysics are thus understanding this relation 
and accurate measurements of cluster mass in observations. 

Dynamical mass estimate uses motions of cluster galaxies 
as tracers of the mass distribution of groups and clusters, and 
it is a relatively inexpensive and efficient technique. The first 
dynamical measu rement of cluster mas s dat es back to the pi- 
oneering work s of lZwickvl d!933l 1 1937b and|Smith| d!936l see 
lBivianoll2000l for a historical overview) . To date, a variety 
of dynamical mass estimation techniques has been devised 
and used extensively, including the projected mass estima- 
tor dHeisler et al.lll985t iRines et all 120031: iRines & Diaferio l 
20061) , the Jeans analysis (ICarlberg et al.lll997t iGirardi et ail 
19981 lyaiid er Mare l et al.l l2000f iBiyiano & Girardil 120031: 
Rines et alj 120031; iKatgert et al.l l2004t iBiviano & Pog gianti 
2009), and the caustic method CDiaferio & Gelled 1 1 997t 
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Rines & Diaferiol 120061: IRines et al l l2007t iDiaferio et aD 
2005ir For low-mass clusters and groups (M < 1O 14 M0), es- 



pecially at higher redshifts, dynamical mass estim ate is often 
the o nly available approach to obtain mass (e.g., ICoil et ail 
2006), since they are mostly X-ray dim and not massive 
enough to be detected by the Sunyaev-Zeldovich effect. Dy- 
namical mass estimates also provide important check on mass 
measurements obtained using other methods an d help in cal- 
ibrati on of optical mass-richness relation (e.g. B ecker et al.l 
12001 . 

Recently, lEvrard et al.l (120081) used a large set of dissipa- 
tionless A^-body cosmological simulations (some with non- 
radiative gasdynamics) and showed that virial scaling relation 
for dark matter particles is consistent with predictions of the 
virial theorem, in which the velocity dispersion scales with 
halo mass to the power of 1/3. Their study further showed 
that the virial scaling relation exhibits remarkable degree of 
universality and self-similarity across a broad range of halo 
masses, redshifts, and cosmological models, as well as hav- 
ing very small scatter. This is good news for the use of optical 
clusters as sensitive cosmological probes. 

In practice, however, applying the virial scaling relation for 
precision cosmology faces a number of difficulties. First, un- 
biased selection of cluster member galaxies (i.e., removal of 
foreground and background interlopers) is difficult and can be 
a major source of systematic uncertainty in ve locity disper- 
sion mesaurement of cluster members (see e. g . , iMcKav et al.l 
120021: IBiviano etal.ll2006t IChen et all 12006). especially for 
low-mass cluste rs and groups consis ting of only a few bright 
galaxies (e.g.. iHernquist et aT1ll995l) . Second, since we can 
only measure the line-of-sight velocity, velocity anisotropy of 
galaxi es is also an additional so urce of systematic uncertainty 
(e.g.. IBiviano & Poggianti 2009). In order to use velocity dis- 
persion as an accurate and precise tracer of cluster mass for 
future optical cluster surveys, it is imperative to understand 
and control these uncertanties. 

However, even if cluster member selection and velocity 
anisotropics are well understood and controlled, we need 
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to understand possible systematic uncertainties in the intrin- 
sic relation between three-dimensional velocity dispersion of 
galaxies and total cl uster mass. For exam ple, the simulations 
used in the study of lEvrard et al.l d2008l) did not include ef- 
fects of dissipative baryonic processes, which do accompany 
cluster formation and lead to formation of cluster galaxies. 
From a number of X-ray observations we know that self- 
similarity is broken in cluster cores and group-size systems 
and this causes deviations of scaling relations from theoreti- 
cal predictions of cluster formation where gravity is the only 
energy driver (e .g., Marke vitchl 1 99 8 ; Arnau d"& Evrardl 19991; 
lOsmond & Ponmanll2004t iPratt et al Il2009l) . Baryon dissipa- 
tion, in particular, modifies the total mass profile and, likely, 
the velocity dispersion profile in groups and clusters both due 
to the contraction an d deepening of cluster potential (e.g., 
Zeldoyich etaLl Il980t iBj umenthaf et al.l 119861; iGnedin et ail 
2004tlSellwood & McGaughll2005l) and due to the ove rall re- 
distribution of baryonic mass (e.g. JRudd et"alll2008l) . Ob- 
servationally, we can only measure velocity dispersions of 
cluster galaxies, which can be a biased tracer of dark mat- 
ter velocity dispersion. This velocity bias is also a poten- 
tial source of systematic uncertainties (e .g., Carlbe rg et alj 
19901: lEvrard"etaiiri994t iFrenk et al.lfl996t Icolfn et al. | 2000f 
Ghigna et al.l20 00; Sprin gel et alj2001l;lpiemand et al. 2004[ 
Gao et al.ll2004lFaltenbacher et alJl2005l:lBiviano et alJbood; 
Faltenbacher & Diemand 2006). 

In this paper, we focus on the latter two intrinsic system- 
atic uncertainties. Namely, we investigate the effect of dissi- 
pation on the dark matter virial scaling relation and velocity 
bias of cluster galaxies using high resolution self-consistent 
cosmological cluster simulations. We show that the baryon 
dissipation increases the velocity dispersion of dark matter, 
by changing the inner density profiles and potential of gravita- 
tionally bound systems, and that the effects are larger in lower 
mass systems where gas cooling and star formation are more 
efficient. We also show that observational estimate of dynam- 
ical mass gives unbiased mass estimate of galaxy groups and 
clusters, in contrary to the previous results obtained based on 
dark matter only simulations. We describe our simulations in 
§[2] and report results in § [3] Main results are summarised in 

§E1 

2. THE SIMULATION 

In this study, we analyze high-resolution cosmological sim- 
ulations of 16 cluster-sized systems performed using the 
Adaptive Refinement Tree (ART) jV-b ody+gasdynamics code 
dKravtsovlll999t iKravtsov et aT][2002h . The ART is a Eule- 
rian code that uses adaptive refinement in space and time , 
and (non-adaptive) refinement in mass dKlvpinet alJ feOOl) 
to reach the high dynamic range required to resolve cores 
of halos formed in self-consistent cosmological simulations. 
The simula t ions pr ese nted here are discus sed in detail in 
iNagai et alj d2007al) and lNagai etail d2007bl) and we refer the 
reader to these papers for more details. Here we highlight 
aspects of the simulations relevant for this work. 

The /V-body+gasdynamics cluster simulations used in this 
analysis include collisionless dynamics of dark matter and 
stars, and gasdynamics. In order to assess the effects of gas 
cooling and star formation on the dynamics of dark matter, 
we conducted each cluster simulation with two different pre- 
scriptions for gasdynamics: one with only the standard gas- 
dynamics for the baryonic component without radiative cool- 
ing and star formation - the 'non-radiative' (NR) runs, and 
'cooling+SF' (CSF) runs. In the CSF runs, several physical 



TABLE 1 

Properties of the simulated CSF clusters at z = 



Cluster ID 


^200c 


lvl 2Q0c 






[ft-l Mpc] 


nn 14 /7 _1 M^i 




CL101 


1 786 


13 2 


1 Z.tJOJt 


CL102 


1 509 


7.98 


749369 


CL103 


1 585 


9.24 


867435 


CL104 


1 435 


6 87 


633740 


CL105 


1 443 


7.00 


649012 


CL106 


1 316 


5.30 


496021 


CL107 


1.174 


3.76 


348465 


CL3 


1.104 


3.13 


981098 


CL5 


0.925 


1.84 


574968 


CL6 


0.958 


2.05 


644496 


CL7 


0.937 


1.91 


603447 


CL9 


0.657 


0.66 


204249 


CLIO 


0.710 


0.83 


261205 


CL11 


0.779 


1.09 


338860 


CL14 


0.755 


1.00 


312626 


CL24 


0.671 


0.70 


221012 



processes critical to various aspects of galaxy formation are 
included: star formation, metal enrichment and thermal feed- 
back due to supernovae Type II and Type la, self-consistent 
advection of metals, metallicity dependent radiative cooling 
and UV hea t ing due to cosmological ionizing background (see 
Nagai et alj d2007al) for details of the metallicity-dependent 
radiative cooling and star formation). These simulations 
therefore follow the formation of galaxy clusters starting from 
the well-defined cosmological initial conditions and capture 
the dynamics and properties of the ICM in a realistic cosmo- 
logical context. However, some potentially relevant physical 
processes, such as AGN bubbles, magnetic field, and cosmic 
rays, are not included. Therefore, the simulated cluster prop- 
erties are probably not fully realistic in the innermost cluster 
regions, where these processes are likely important. 

Our simulated sample includes 16 clusters at z = and their 
most massive progenitors at z = 0.6 and z= 1.0. The proper- 
ties of simulated clusters at z = are given in Table Q] The 
total cluster masses are reported at the radius r2ooc enclosing 
overdensities with respect to the critical density at the red- 
shift of the output (below, we also use a higher overdensity 
level, 500). The dark matter (DM) particle mass in the re- 
gion around the cluster was 9 x 10 8 /T 1 M Q for CL101-107 
and 3 x 10 8 /j _1 Mq for CL3-24, while other regions were 
simulated with lower mass resolution. The peak resolution 
of the cluster simulation is w 2 - 3 h kpc in the box size 
of 80-120/T 1 Mpc. The simulations are performed for the 
flat ACDM model: fl m = 1 - fl A = 0.3, Q b = 0.04286, h = 
0.7 and a& = 0.9, where the Hubble constant is defined as 
100/; km s" 1 Mpc" 1 , and an erg is the power spectrum normal- 
ization on an 8/i _1 Mpc scale. 

3. RESULTS 

3.1. Virial Scaling Relation 

Following lEvrard et alJ (12008). we assume velocity 
isotropy and calculate the one-dimensional dark matter veloc- 
ity dispersion cr DM to be 

j 3 N 

;=i p=i 
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FIG. 1 . — Left panel shows the dark matter virial scaling relation for the sixteen clusters for A c = 200 using least-square fitting at z = and z = 1 and the right 
panel shows the same relation for A c = 500. For each panel, the CSF clusters at z = are represented by the solid blue circles, and the CSF clusters at z = 1 are 
represented by open blue circles. The solid blue line is the best least-square fit for the z, = 0. For the NR ran, the z = and z = 1 clusters are represented by solid 
red squares and op en red squares respectively, which the red dot-dashed line is the best least-square fit to the clusters at z, = 0. For reference, the best fit from 
Evrard et al. (2008) is shown as black dashed line. Please refer to Table|2]for the values of the fitted parameters. 

highest binding energy. 

A power law is then fitted to the dark matter virial scaling 
relation using standard least square minimization: 

Kz)M Ac 



TABLE 2 
Scaling Relation Parameters 







A c = 200 








z = 0.0 


z = 0.6 


Z= 1.0 


a 




NR 


0.3072 ±0.01 14 


0.3214 ±0.0246 


0.3512±0.0152 






CSF 


0.2724 ±0.0149 


0.2903 ±0.0254 


0.3232 ±0.0153 


°T>M,0 




NR 


643 ±7 


648 ± 1 1 


645 ±8 


[km/s] 




CSF 


692 ± 1 1 


674 ± 1 1 


671 ±8 


<*L> 


n 


NR 


0.0410 ±0.0079 


0.0630 ± 0.0086 


0.0348 ± 0.0062 






CSF 


0.0557 ±0.0092 


0.0569 ±0.0098 


0.0375 ± 0.0080 






Ac = 500 








z = 0.0 


Z = 0.6 


z= 1.0 


a 




NR 


0.3335 ±0.0096 


0.3256 ±0.0207 


0.3475 ±0.01 10 






CSF 


0.2965 ±0.0144 


0.2888 ±0.0180 


0.3134±0.0108 


°TJM.O 




NR 


742 ±6 


737 ± 10 


751 ±10 


[km/s] 




CSF 


788 ±10 


762 ±9 


771 ±11 




/i 


NR 


0.0310 ±0.0059 


0.0444 ± 0.0082 


0.0258 ± 0.0062 




CSF 


0.0459 ±0.0075 


0.0394 ± 0.0073 


0.0361 ±0.0070 



log <r DM (M, z) = log cr DM ,o + a log 



2 X 10 14 Mf; 



(2) 



where N is the total number of dark matter particles inside 
r& c , the radius which encloses matter having average density 
A c times the critical density; v ; ,., is the the i-th velocity com- 
ponent (i = x,y,z) of the p-th dark matter particle inside r& c ; 
and v, is the peculiar velocity of the cluster is defined to be 
the mass-weighted average of all dark matter particles within 
?A C - We have also calculated the cluster velocity using veloc- 
ities of all species (dark matter, stars and gas), and it agrees 
with the average velocities of dark matter to within a few per- 
cent for most clusters, since baryons are only a minor fraction 
of the total cluster mass. We weigh the mean and variance by 
dark matter particle mass but it is insensitive to whether mass- 
weighted or number- weighted is used within r&, c . We take the 
cluster center to be the location of the subhalo which has the 



where h(z) = H(z)/100 km s" 1 Mpc" 1 is the normalized Hubble 
parameter which takes into the account of the redshift depen- 
dence of Ma c . Figure Q] shows the fitted velocity dispersion- 
mass relations of the CSF and the NR runs at z = 0, 0.6, and 
1 .0. Table [2] reports the best fit parameters of the velocity 
dispersion-mass relations measured at A c = 200 and 500. The 
errors are estimated by generating 5000 bootstrap samples 
from the sixteen clusters at each redshift. 

At A c = 500, the NR scaling relation is consistent with the 
self-similar scaling relation with no evolution. At A c = 200, 
the best-fitted NR scaling relation at z = differs slightly 
(2.6cr) from that of lEvrard et ail (12008b with a = 0.3361 ± 
0.0026 and cr D M(M 20 o c = lO^/z^Mg) = 1082. 9 ±4.0 km/s 
Note t hat results are fully consistent with the lEvrard et alj 
(2008) relation if we limit our analyses to systems with the 
mass range (M > 10 14 Ii~ x Mq) adopted in their analysis. 
Larger number of simulated clusters is needed to explore re- 
maining minute deviation from the self-similar relation. 

The CSF scaling relations differ significantly from the NR 
relation at all redshifts. The baryon dissipation increases the 
velocity dispersion of dark matter, making the slope shallower 
in the CSF run compared to that of the NR run. Interestingly, 
the log-normal scatter is only 3-5%, and it is remarkably in- 
sensitive to input gas physics and redshift. 

The breaking of self-similarity by baryon dissipation in the 
CSF runs is illustrated more clearly in Figure 12 which shows 
the ratio of dark matter velocity dispersion in the CSF run 
relative to the NR run, o^Mc / a 2W)c 5 as a function of halo 

5 As the same cluster in CSF and NR runs may have slightly differ- 
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FIG. 2. — The ratios of velocity dispersions in the CSF and NR runs, 
CT 200c/ <J 200<: as a f uncnon °f cluster mass at z = 0. Different symbols indi- 
cate velocity dispersions computed using DM particles in different annulus: 
[0-1] (circles), [0.2-1] (square) and [0.4-1] (triangles) in unit of r2ooc- 

mass at z = on a cluster-by-cluster basis. For almost all 
halos, the inclusion of baryonic physics leads to an increase 
in dark matter velocity dispersion, and the effect is greater 
for lower mass systems. We find that this bias is of the or- 
der ~ 10%- 15% for M 20 o c ^6x IO'VMq, compared to 
6% for M200C > 3 x 10 14 !i~ 1 Mq. This is because gas in lower 
mass halos is at a lower virial temperature and thus capable 
of cooling more efficiently to form stars. We find the same 
trend at higher redshift. Note that there is one outlier (CL5) 
which shows a decrease in the velocity dispersion in the CSF 
run. This is due to fast moving subhalos in the NR run that 
have already merged with the central subhalo in the CSF run. 
Another outlier in the opposite direction (CL6) has a ~ 20% 
increase in velocity dispersion. This is also due to a similar 
merging event in the CSF, boosting its velocity dispersion. 

We further show that the change in velocity dispersion is 
largely due to the baryon dissipation in the central regions 
of clusters. To illustrate this, we also compute velocity dis- 
persions by excluding the central regions of clusters in Fig- 
ure |2] This shows that differences in velocity dispersions in 
the CSF and NR runs become smaller when we apply central 
exclusion. Excluding r < 0.2r2oo reduces the differences in 
the CSF and NR to about 4% on most systems. The bias is 
even smaller (less than 2%) on all systems if we exclude the 
central r < 0.4r2oo- This illustrates that the increase in <7200c is 
mainly driven by the increase in velocities in the inner regions 
of CSF groups and clusters. 

Figure [3] shows the velocity dispersion profiles averaged 
over the 16 clusters at z = and z = 1 for both CSF and NR 
runs. Note that the mass depedence seen in Figure [2] is ex- 
plicity taken out by dividing the circular velocity v 2 00c The 
dark matter velocity dispersion of CSF clusters increases to- 
wards smaller radii compared to their NR counterparts. The 
change in the shape of the velocity dispersion profile is a di- 

ent M2000 we normalize each a^ooc w ' m tne crrcular velocity V200C = 
\J GM2ooc/r200c- 



0.1 1 0.1 1 

r/r 200c r A 2 ooc 

FIG. 3. — Average DM velocity dispersion profiles for the sixteen clusters 
at z = (left-panels) and z = 1 (right-panels) for CSF (blue solid lines) and NR 
(red dashed lines) runs. Dotted lines are the corresponding la errors. The top 
panels show 1-D differential velocity dispersion profiles; the middle panels 
show the 1-D cumulative velocity dispersion profiles; the bottom panels show 
the corresponding number density profiles. Also plotted are profiles for CSF 
cluster galaxies (blue star), CSF subhalos (blue circle) and NR subhalos (red 
squares). 

rect consequence of the enhancement of dark matter density 
caused by adiab atic contraction repo rted previously based on 
our simulations dGnedin et al.l 120041) : as gas cools and con- 
denses towards the center, it deepens the central potential and 
drags dark matter inward. 

Note that we have assumed dark matter velocity dis- 
perion to be isotropic. This, however, is not a com- 
pletely valid assumption. In fact, some level of velocity 
anisotropy is expected for collisionless dark matter particles 
dCole & Lacevlfl996l iCarlberg et al J [19971: iCohn et aUl2000t 
lHansen & Moordl2006l) . Cluster galaxies are also collision- 
less as they have relaxation time greater than the Hubble time. 
Velocity anisotropy is therefore an important ingredient in 
the dynamical mass estimates. The anisotropy is commonly 
quantified by the anisotropy parameter /3, defined as 

P= l-OSaf/aj, (3) 

where 07 and 07 are the radial and tangential velocity disper- 
sions respectively. The definition of a follows from Eq. ([TJ), 
where contribution of mean velocities, which is non-zero be- 
yond the virial radius, is included. 

We therefore use our simulations to study the velocity 
anisotropy of dark matter as well as galaxies in groups and 
clusters. Figure [4] shows that dark matter velocity is nearly 
isotropic in the inner region of halo at 0. 1 r200c, and it becomes 
increasingly radial at larger radii (e.g., j3 = 0.3 at 0.5r2ooc and 
= 0.35 at r2ooc)- Interestingly, we find that the dark matter 
velocity anisotropy profile is almost unaffected by the addi- 
tion of cooling, star formation and feedback and insensitive to 
redshift between z = and 1 . The CSF velocity anisotropy is 
only slightly lower than that of the NR one. This makes sense 
given that the velocity anisotropy profile is o nly dependent on 
the sl ope of the density profile at large radii (lHansen & Moord 
2006), where the effect of baryon dissipation is small. 
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FIG. 4. — The average velocity anisotropy profiles. The notations are the 
same as FigurefS] 



Finally, lEvrard et alJ ([2008) show that low resolution can 
bias the slope of the scaling relation high. Their resolution 
studies indicate that the calibration of the virial scaling rela- 
tion is not affected by resolution if clusters are resolved with 
more than 10 5 dark matter particles. In this work, we ensured 
that all groups and clusters used in our analyses are resolved 
with greater than 3 x 10 5 particles. Thus, the resolution of our 
simulations is sufficient for the virial scaling relation work. 
The magnitude of the baryon dissipation, on the other hand, 
depends sensitively on the cooled gas fractions in simulations. 
Given that the current simulations suffer from the overcooling 
problem (e.g jBorgani & Kravtsovl20 09). the effect of baryon 
dissipation in this work should be taken as an upper limit. 

3.2. Velocity Bias 

Since dark matter velocities are not directly observ- 
able, galaxies are used as tracers of the dynamics of 
underlying dark matter. Velocity dispersion of clus- 
ter galaxy er ga i, however, is believed to be a "biased" 
tracer of the dark matter dynamics, opm dCarlberg et al] 
19901 lEvrard et al.l|1994tlFrenk et alJll996tlColfn et alJl2000h 
Ghigna et al.1 l2000l: ISpringel et al] fcOqiTlR dtenbache r et al] 



Faltenbacher & Diemand 2006; Benatov et al. 2006) 



This velocity bias is often defined as, 

0"DM 



(4) 



where, as before, the one-dimensional velocity dispersions 
of both galaxies and dark matter are calculated assuming 
isotropy (i.e., Eq. (fl}). 

Figure [5] show s the velocity bias f or bo th the CSF and NR 
runs. Following iNagai & Kravtsovl d2005l) . we select subha- 
los using their total bound mass (M tot > 3 x 10 10 h~ l M Q ) and 
stellar mass (M» > lO'/r'M©). In the NR runs, we only have 
the first sample, since there are no stars in these simulations. 
Results are summarized in the Table [3] 

A main finding is that the velocity bias depends on how 
galaxies are selected. In the NR runs, we select subhalo us- 
ing their current total bound mass and find that there is a 
positive bias of about 10% at all redshifts. This is in good 
agreement with the positive velocity bias of w 10% - 15% 
for dark matter subhalos in jV-body simulati ons dColfn et al] 
I2000t iDiemand et all l2004t lGao et al.1 [2004). This positive 
velocity bias arises because selection of subhalos using their 
current bou nd mass radially biases the resulting sample to 
larger radii (Nagai & Kravtsovll2005l) . This selection also bi- 
ases the sample by including more recently accreted subha- 
los which have larger velocity dispersions. Tidal disruption 
can lead to not only mass loss but also complete disruption 




0.1 1 
h(z)M 800c [10"»M o ] 



0.1 1 
Mz)M 500c [10«>M o ] 



FIG. 5. — Left panel shows the velocity bias for A c = 200 at z = 0, while 
the right panel shows the same for A c = 500. For each panel, the red squares 
represent the velocity bias for the NR subhalos selected using their current 
total bound mass of > 3 X 10 /t'Mq, the blue circles are velocity bias for 
CSF subhalos selected with the same mass threshold. Blue stars show the 
velocity bias for the CSF subhalos selected using their stellar masses with 
threshold M* > 10 9 Ii~ 1 Mq. We refer reader to the Table [5] for the mean 
velocity bias and 1 — a scatter at higher redshifts. 



(or, equivalently, to mass loss that brings subhalo mass be- 
low resolution of simulation). The disrupted subhalos tend to 
be biased towards smaller radii and cannot be recovered by 
any selection criterion. It is generally thought that dissipa- 
tion should enhance the survival of subhalos, although recent 
simulations sho w that effect is not as large as has been previ - 
ously thought dNagai & Kravtsovl 120051 iMaccio et al.ll2006t) . 
Indeed, the numbers in Table[3]show that using the same sub- 
halo selection criteria in the CSF runs the velocity bias is no- 
ticeably smaller (by a few percents). This means that baryon 
dissipation makes subhalos somewhat less susceptible to tidal 
disruption in the inner regions of groups and clusters. 

In the CSF runs, we also select galaxies based on their stel- 
lar mass. Interestingly, galaxies selected this way shows al- 
most no velocity bias. The selection criteria are important 
here because stars formed in the inner regions of galactic-size 
halos are tightly bo und to the system and ar e more resistant 
to tidal mass loss (Nagai & Kravtsov 2005). This also ex- 
plains why the bias is generally smaller at higher redshifts, 
where galaxies has not yet h ad the time to exper ience sig- 
nificant tidal mass loss (e.g., IConroy et al.l 120061) . This is 
also consistent with the previous results that find that both 
radial and velocity biases of subhalos largely disappear if 
subhalos are selected using their p roperties at the accretion 
epoch, unaffected by tidal strip ping (Nagai & Kravtsovll2005l : 
iFaltenbacher & Diemand 2006). Moreover, this explains why 
the bias is slightly higher in the inner region (e.g., A c = 500), 
where the tidal effects should be more significant. 

The importance of how subhalos are selected is further 
highlighted in Figure [3] which compares galaxy velocity dis- 
persion profiles (points with errorbars) to the velocity disper- 
sion profiles of dark matter particles, indicated by lines. The 
figure shows that velocity dispersions of galaxies are gen- 
erally higher than those of dark matter in the radial range 
0.1 < r/r2ooc < 1- The velocity bias is slightly larger when 
galaxies are selected using their current total mass, rather than 
their stellar mass. The figure also illustrates that both velocity 
and spatial biases are smaller at higher redshifts, where the 
tidal stripping and disruption had not yet had time to affect 
the masses of subhalos, and hence their spatial distribution 
and overall velocity dispersion. 

For high-redshift objects, the number of cluster galaxies 
Afgai is often too limited for reliable dynamical mass. Using 
CSF simulations, we quantify how well one can estimate dark 
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TABLE 3 
Velocity Bias 



mean±l-<r error 


A c = 200 






z = 0.0 


z = 0.6 


z.= 1.0 






1.131 ±0.022 


1.063 ±0.036 


1.073 ±0.032 






1.067 ±0.021 


1.050 ±0.021 


1.032 ±0.024 


fe dSF 




1.029 ±0.022 


0.993 ±0.032 


0.971 ±0.033 




3) 


0.796 ±0.070 


0.792 ±0.055 


0.772 ±0.065 




5) 


0.857 ±0.042 


0.783 ±0.053 


0.792 ± 0.047 




10) 


0.964 ±0.035 


0.894 ±0.027 


0.856 ±0.047 




20) 


0.976 ±0.025 


0.969 ±0.027 


0.957 ±0.037 






A c = 500 






z = 0.0 


z = 0.6 


z= 1.0 


^\ 

V.Sltl) 




1.200 ±0.021 


1.084 ±0.043 


1.1 16 ±0.039 


v.sub 




1.1 28 ±0.028 


1.084 ±0.021 


1.078 ±0.025 


b CSF 

v,gal 




1.083 ±0.031 


1.018 ±0.034 


1.020 ±0.028 


*SS ^ga. = 


3) 


0.811 ±0.077 


0.756 ±0.066 


0.816±0.067 


«(^ = 


5) 


0.945 ± 0.052 


0.821 ±0.052 


0.887 ±0.057 


*S (iV ^ = 


10) 


1.013±0.043 


0.928 ± 0.037 


0.958 ±0.039 


= 


20) 


1.030 ±0.035 


1.022±0.019 


1.067 ±0.037 



matter velocity dispersion using a limited number of galaxies. 
Specifically, we compute the velocity bias using 3, 5, 10 and 
20 most massive galaxies (excluding the central galaxy) in our 
CSF halos. The results are summarized in Table[3] Our analy- 
sis shows that galaxy velocity dispersion tends to exhibit less 
bias for larger number of galaxies. At z = 0, galaxy velocity 
dispersion is essentially unbiased (i.e., b v = 1) within lex for 
yVg a i > 10 at both A c = 200 and 500. At higher redshifts, larger 
number of galaxies are required to achieve unbiased estimate. 
Our analyses indicate that the galaxy velocity dispersions are 
unbiased within 2<r for YV ga i > 20 at z = 0.6 and 1. A larger 
sample of simulated clusters are required to confirm the re- 
maining bias at a level of 5% at these redshifts. Moreover, 
our results clearly show that the velocity bias is significant for 
A^ a i < 10 galaxies. This is because the more massive galaxies 
reside nearer to the inner region due to dynamical friction and 
hence have lower velocities. 

4. SUMMARY AND DISCUSSION 

In this work we investigate the effect of baryon dissipation 
on the scaling relation between total mass and velocity dis- 
persion, and the velocity bias of galaxies in clusters. To this 
end, we use a representative set of self-consistent cosmologi- 
cal cluster simulations that follow the formation and evolution 
of galaxies by taking into account effects of gas cooling, star 
formation, and feedback. 

We confirm that the virial scaling relation of dark matter ex- 
hibits a remarkable degree of self-similarity with log-normal 



scatter of only 5%. The virial scaling relation in our non- 
radiative hydrodynamical simulations is consiste nt with that 
of the dissipationless simulations presented by lEvrard et al.l 
(2008). This demonstrates that the virial scaling relation has 
now been calibrated at a level of a few percent in the non- 
radiative regime. 

We show that the largest source of systematic uncertainty in 
the intrinsic virial scaling relation calibration lies in our un- 
derstanding of baryonic physics. Using numerical simulation 
that include the physics of galaxy formation, we find that the 
baryon dissipation increases dark matter velocity dispersion 
of of groups and clusters at the level of ~ 4% to 10% due 
to adiabatic contraction of dark matter in the central regions. 
Our simulations indicate that this effect is more pronounced in 
lower mass systems where gas cooling and star formation are 
more efficient than massive systems. Note, however, that the 
effect of baryon dissipation in this work should be taken as an 
upper limit, because the c urrent simulations suffer fro m the 
overcooling problem (e.g.. iBorgani & Kravtsovl 12009). De- 
spite this problem, we find that the scatter in the virial scaling 
relation is very tight (only ~ 5%), independent of mass and 
redshift, and insensitive to details to input cluster physics. 

In practice, the dark matter virial scaling relation is not 
direclty accessible to observations. Instead, observers rely 
on galaxies as tracers of the underlying distribution and dy- 
namics of matter. Thus, the bias in the velocity dispersions 
of galaxies with respect to that of dark matter is one of the 
main sources of systematic uncertainty in the mass estimates 
of galaxy clusters. Previous works based on TV-body sim- 
ulations have shown that dark matter "subhalos" in simula- 
tions exhibits a positive velocity bias, because slow subhalos 
are much less common, due to physical disrupti on by grav- 
itational tides early in the merging history (e.g., IColrn et all 
l2000llDiemand et al.ll200l: iGao i et al.ll2004h . 

In this work, we show that the galaxy velocity bias, in fact, 
depends on how galaxies ("subhalos") are selected in simu- 
lations. Using the self-consistent cluster simulations, we find 
that galaxies selected based on their stellar mass exhibit no 
velocity bias. Moreover, we find that the observational es- 
timates of galaxy velocity dispersion with more than twenty 
most massive galaxies in groups and clusters is an unbiased 
tracer of the dynamics of underlying matter. 
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